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Abstract 

Cosmic-ray photons above 10^^ eV can convert in the geomagnetic field and initi- 
ate a preshower, i.e. a particle cascade before entering the atmosphere. We compare 
the preshower characteristics at the southern and northern sites of the Pierre Auger 
Observatory. In addition to a shift of the preshower patterns on the sky due to 
the different pointing of the local magnetic field vectors, the fact that the northern 
Auger site is closer to the geomagnetic pole results in a different energy dependence 
of the preshower effect: photon conversion can start at smaller energies, but large 
conversion probabilities (> 90%) are reached for the whole sky at higher energies 
compared to the southern Auger site. We show how the complementary preshower 
features at the two sites can be used to search for ultra-high energy photons among 
cosmic rays. In particular, the different preshower characteristics at the northern 
Auger site may provide an elegant and unambiguous confirmation if a photon signal 
is detected at the southern site. 



1 Introduction 



Substantial fluxes of cosmic-ray pliotons at ultra- high energy (UHE), above 
10^^ eV, are predicted by non-acceleration (top-down) models of cosmic-ray 
origin (for example, see [1]). At a smaller level, UHE photons are also expected 
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to be produced in acceleration (bottom-up) models [2]. UHE photons can be 
detected through the particle cascades they initiate in the atmosphere. So 
far, upper limits on the photon flux were set PP1I51I51I7] . The large exposure 
expected to be collected during the next years, in particular by the Pierre 
Auger Observatory [8J, will enormously increase the sensitivity for detecting 
UHE photons [7J. 

Contrary to the case of hadron primaries, UHE photons around 10^'^ eV can in- 
teract with the magnetic fleld of the Earth before entering the atmosphere [9J. 
This process is commonly called geomagnetic cascading or preshower and 
leads to a dramatic change of the air shower development for primary pho- 
tons pTOllTTllTMMi] . 

The probability of magnetic e^e~ pair production ( "photon conversion" ) and, 
in case of conversion, the synchrotron emission by the produced electrons de- 
pend on the particle energy and on the transverse component of the local 
magnetic fleld p^Sfll] . Photon conversion can occur at altitudes up to sev- 
eral thousand kilometers above the ground, i.e. well above the atmosphere 
(~100 km altitude). A preshower consists mostly of the secondary photons 
emitted by the electrons. A secondary photon of sufficiently high energy may 
additionally convert in the magnetic fleld and produce another electron pair. 
The number of preshower particles depends on the speciflc conditions and can 
fluctuate strongly; it can range from a few particles to a few thousand. Nev- 
ertheless, as we demonstrate in this paper, the preshower formation can be 
fairly well inferred from characteristics of the observed extensive air shower. 
The spread of the preshower particles in transverse distance from the shower 
axis and in arrival time at the atmosphere is well below the resolution of cur- 
rent shower experiments [H]. Hence, also in case of geomagnetic cascading, 
the subsequent air shower would be observed as one event by the experiments. 

Since the probability of photon conversion and the synchrotron emission spec- 
trum depend on the local transverse magnetic fleld, characteristics of the 
preshower process (such as conversion probability, or the energy spectrum of 
the preshower particles) depend strongly on the speciflc trajectory through the 
magnetosphere. For a given experimental site, this leads to a large directional 
dependence of preshower characteristics within the local coordinate system of 
zenith 9 and azimuth 0. In addition, as already noted in the work of McBreen 
and Lambert 25 years ago [9j, a dependence is expected on the location of the 
experiment: experimental sites at different geomagnetic latitudes have differ- 
ent local magnetic fleld conditions, which can affect preshower characteristics 
for UHE photons. 

The comparison of preshower characteristics of UHE photons expected at two 
sites located at different hemispheres of the Earth is of particular current inter- 
est due to the construction of the Pierre Auger Observatory. The southern part 
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of the Observatory ( "Auger South" ) is situated near Malargue (Argentina) at 
69.2° W, 35.2° S. The construction of Auger South is approaching completion 
and first results were obtained, including an upper limit to UHE photons [7J. 
To achieve full sky coverage, the northern site ( "Auger North" ) is planned in 
Colorado (USA) at 102.7° W, 37.7° N. At Auger South, the magnetic field of 
~24.6 fiT points upward to ~ 55°, ~ 87°. At Auger North, the magnetic 
field of ~52.5 jiT points downward from 6 ~ 25°, (j) ~ 262°{3 

Previous works in the literature focused rather on the directional dependence 
of preshower characteristics of UHE photons for one given site than compar- 
ing the conditions at different sites. Examples of magnetic field strengths at 
various sites were given in [11]. In that work, the observable of shower size 
(total number of electrons at ground) was studied. It was noted that azimuth 
regions of smaller (or stronger) preshower effect can be identified that may 
differ between the sites. In [12], the altitude dependence of the transverse 
component of the magnetic field was compared for the two Auger sites for 
a few directions, and preshower particle spectra were studied. E.g. in [13], it 
was pointed out that at Auger North, due to the larger magnetic field, the 
preshower effect starts at a lower energy. A particular focus of that work was 
the study of observables for the ground array, where differences in the detector 
signals for converted and unconverted primary photons were found. 

In this work, making use of our PRESHOWER code p3] linked to CONEX, 
a fast shower simulation program [16], we perform a dedicated comparison 
of the preshower characteristics of UHE photons for the conditions of the 
Auger North and Auger South sites. The geomagnetic field in PRESHOWER 
is described by the IGRF model [T7IITB] . We find that the transition range in 
energy from small to large conversion probabilities is larger at Auger North, 
starting at smaller but ending at higher energy. We study how the differences 
of the preshower characteristics between the two sites can be used to perform 
a complementary search for UHE photons. In particular, the situation may 
arise where a possible UHE photon signal detected at Auger South requires an 
independent cross-check. As will be shown, such a cross-check could come from 
Auger North. In fact, making use of the UHE preshower effect by observing 
from a different experimental site seems the most elegant and natural way for 
cross-checking an UHE photon signal claim. 

While the study is performed for the specific case of the two Auger sites, the 
general findings hold for any two sites with sufficiently different local magnetic 
field conditions. The results shown for Auger North hold, to a good extend, 
also for the HiRes [19] and Telescope Array [20] sites, both of them being 
located in Utah. 



Azimuth is defined in this work counterclockwise from geographic East. For in- 
stance, = 0° means East, (j) = 90° North etc. 
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The plan of the paper is as follows. In Section [21 preshower features are dis- 
cussed, and in Section [3l the corresponding air shower features are studied. 
In Section m UHE photon scenarios and their observation at the two different 
sites are investigated. In particular, it is worked out how observations at the 
two sites can complement each other. Conclusions are given in Section [51 In 
Appendix [X]it is shown how the transverse component of the geomagnetic field 
in the IGRF model varies along different particle trajectories. In Appendix [Bl 
we show a collection of conversion probability sky maps supplementing Section 

urn 



2 Comparison of preshower features at Auger South and North 

2.1 Conversion probability 

A key parameter to characterize the fate of an UHE photon in the Earth's 
magnetic field is the conversion probability Pconv Given the local differential 
probability of a photon to convert into an electron pair, Pconv results from an 
integration along the particle trajectory. Small values of Pconv indicate a large 
probability of the UHE photon to enter the atmosphere without conversion 
and to keep its original identity. In turn, UHE photons would almost surely 
undergo geomagnetic cascading for values of Pconv close to unity. 

Pconv depends on the experimental site, the photon energy, and the direction of 
the particle trajectory in the local coordinate system of zenith 6 and azimuth 
0, Pconv = f{s^iQ,E,0,(f)). Thus, for a chosen site and a fixed primary photon 
energy, sky maps within the local coordinate system Pconv = fi^, 0) can be 
produced to study the pattern of UHE photon conversion. 

As an example, in Figure [H sky maps are shown for two different energies 
lg(E/eV) = 19.6 and 20.0, corresponding to ~ 40 EeV and 100 EeV, both for 
Auger North and South (sky maps at various energies between lg(i?/eV) = 19.4 
and 20.4 are collected in Appendix [B]). One can see that at 40 EeV, photon 
conversion is already important at Auger North for a large part of the sky. 
In contrast to this, at Auger South conversion is very unlikely for 40 EeV 
photons independent of their arrival direction. 

The energy above which photon conversion becomes important, is expected 
to be smaller at Auger North due to the larger local magnetic field strength. 
To quantify this, we show in Figure [2l a sky map of photon energies for a fixed 
conversion probability of Pconv = 10%. It can be seen that photon conversion 
starts at a factor ~2 higher energy at Auger South compared to Auger North. 
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Fig. 1. Sky maps of conversion probabilities at lg(ii^/eV)=19.6 (left) and 20.0 
(right) for Auger North (AN, top), Auger South (AS, middle), and the difference 
(AN— AS, bottom). Contour lines are given for the conversion probability (top and 
middle rows of plots) with a stepsize of 10%, and for the difference of conversion 
probabilities (bottom plots) with a stepsize of 20%. Azimuthal directions are labeled 
("E" for East etc.). Zenith angles are given as concentric circles of 10° steps {9 = 0° 
in the center). The pointing direction of the local magnetic field vector at ground 
is indicated for a specific site (top and middle plots). 



At energies of 100 EeV, also for Auger South photon conversion is important 
for a sizeable part of the sky (Figure [1], right). As expected, small conversion 
probabilities are found for sky regions around the pointing direction of the 
local magnetic field vector. Thus, we expect significant differences in Pconv 
between the two sites for a given direction in terms of local coordinates, with 
the difference being strongly dependent on the regions of the sky. Also shown 
in Figure [H (bottom plots) are therefore sky maps of the difference in Pconv 
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between Auger North and South for the two fixed energies. While for the larger 
part of the sky, Pconv is larger at Auger North (due to the stronger field), there 
is at 100 EeV a significant part of the sky where Pconv is larger at Auger South 
(due to the different pointing of the magnetic field vector). 

It is clear from this that cuts on the local shower arrival direction can be 
introduced to select regions of the sky where Pconv is larger (or smaller) at 
one site compared to the other site. Applying the same cuts to data from the 
two sites, event sets with enriched (for the one site) and depleted (for the 
other site) fractions of converted (or unconverted) photons can be prepared. 
A possible photon signal could then show up with different signatures at the 
two sites for the same selection cuts. This will be studied in more detail in 
Section HI 

Another feature is showing up in Figure [2] where sky maps of energies are 
plotted for fixed Pconv The range of energies covered (for the same value of 
Pconv) is larger at Auger North than at Auger South. For instance, the energy 
at which Pconv = 10% is reached changes at Auger North by a factor ~15 over 
the sky, at Auger South only by a factor ~4. In other words, the sky pattern 
for geomagnetic cascading is more "inhomogeneous" at Auger North (larger 
transition range in energy from small to large conversion probability) than at 
Auger South. 

In fact, while the preshower effect starts at smaller energies for a larger part of 
the sky at Auger North, the transition range to full conversion extends to even 
higher energies at Auger North, see Figure [2] (bottom), sky map for Pconv = 
90%. As an example, a 250 EeV photon (lg(i?/eV)=20.4) has a conversion 
probability > 40% for all directions at Auger South, while at Auger North a 
small window with Pconv < 10% still exists (see Figure [R3l) . 

At first glance this seems unexpected. As the magnetic field strength is a fac- 
tor two larger at Auger North, one would rather expect photon conversion 
at smaller energies at Auger North. However, this feature is again related to 
the different locations within the geomagnetic field, with Auger North being 
considerably closer to the magnetic pole: in this case, the local magnetic field 
lines are less curved when following them from ground to higher altitudes. 
Correspondingly, the transverse projection of the magnetic field vector, which 
determines the local conversion probability, remains small out to considerably 
larger distances from the Earth compared to the situation at Auger South. 
Thus, the fact that Auger North is closer to the pole results in a larger tran- 
sition region for photon conversion, starting at smaller energies compared to 
Auger South (due to the larger strength of the local magnetic field) and ending 
at higher energies compared to Auger South (due to the smaller bending of 
the magnetic field lines with distance). 
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AN, P=10% 




AS, P=10% 




Fig. 2. Sky maps of photon energies for fixed conversion probabilities of 10% (top), 
50% (middle) and 90% (bottom) for Auger North (left) and Auger South (right). 
Contour lines are given with a stepsize of Alg(£^/eV) = 0.1. Values of the mini- 
mum and maximum energies are assigned to the lines (in lg(£^/eV)). See caption of 
Figure [1] for further explanations. 



To summarize, the different positions of tlie sites in tlie geomagnetic field 
translate into a spatial shift of the conversion patterns at the sky according 
to the pointing of the local magnetic field vector and into modifications of the 
transition energy range. The modifications of the transition energies are not 
a simple shift but also the size of the energy range is affected. 
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Fig. 3. Average energy spectrum of preshower particles, weighted by particle energy 
(converted events only) for 10^*^ eV photons. Chosen are two directions, one at each 
site, with the same integrated conversion probability Pconv = 50%. Direction 1 
(Auger North): ~ 60°, </) ~ 270°, direction 2 (Auger South): ~ 66°, ~ 150°. 



2.2 Particle energy spectra in converted showers 



Here we check to what extent it is sufficient to consider only the difference 
in conversion probabihty when comparing the two sites. For a given energy, 
two different directions with the same integrated conversion probability Pconv 
correspond to different spatial trajectories through the magnetosphere. One 
could imagine in one case a relatively small transverse magnetic field extend- 
ing over a larger distance, while in the other case relatively close to the Earth, 
a stronger magnetic field contributes. The same integrated conversion prob- 
ability could be reached, but the spectrum of preshower particles may differ 
(for instance, the synchrotron emission spectrum gets harder with stronger 
magnetic field). 

In Figure [3|, averaged spectra of preshower particles at the top of the atmo- 
sphere are shown, weighted by the particle energy, for different directions at 
Auger North and South. In both cases, Pconv = 50%. One can see that both 
spectra look quite similar; most energy is carried by particles of ~ 10^^ eV. 
A shift between the distribution of about a factor 1.6 in energy can be noted, 
however. Performing air shower simulations for these two directions, a differ- 
ence in the average depth of shower maximum X^ax of the converted events 
of ~ 20 g cm~^ is obtained. This difference is much below the shift in X^ax 
comparing converted and unconverted photon showers to each other (see next 
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Section). For other directions of same integrated conversion probability, dif- 
ferences are typically smaller than shown here. For our purposes in this work, 
we conclude that the differences in conversion probabilities are a sufficient 
measure for the different preshower characteristics when comparing the two 
sites. 



3 Air showers initiated by converted and unconverted photons 

In this section, we study the basic features of air showers initiated by converted 
and unconverted photons, with particular emphasis on the differences between 
these two classes of events. The reasoning is that the differences between Auger 
South and North in the sky patterns for photon cascading, studied in the 
preceeding Section, translate into differences of the air shower features that 
are well observable. 

It is well known that unconverted showers, contrary to converted, have a con- 
siderably delayed shower development due to the LPM effect [H]. Addition- 
ally, event-by-event fluctuations can be extraordindarily large due to a positive 
correlation of the suppression of the cross-section with air density (a photon 
traversing a certain depth AX without an interaction has a larger probability 
for traversing an additional AX due to the larger air density). 

In Figure m examples of shower profiles are shown for 10^'' eV primaries. Sim- 
ulations are carried out with CONEX, which reproduces well CORSIKA [22] 
results. The extrapolation of the photonuclear cross section is based on the 
fit by the Particle Data Group [2B]. All primaries arrive from the same lo- 
cal direction, in this case from geographic North at 45° zenith. The photon 
conversion probability is large at Auger North (> 99.9%) and small at Auger 
South (0.4%). Consequently, most of the photons at Auger North convert and 
have a depth of shower maximum 200—300 g cm~^ smaller than the (mostly 
unconverted) photon showers at Auger South. As expected, also fluctuations 
are much smaller at Auger North in this example. The corresponding distri- 
bution of depth of shower maximum Amax from simulating 1000 events for the 
conditions in Figure H] is shown in Figure [5] (upper panel). For the same local 
direction, the expected features of photon showers are very different at the 
two sites. 

For comparison, also profiles for hadron primaries are added using the QGS- 
JET 01 model |23]. Converted photons are more similar to protons, but still 
reaching A^ax at larger depths. More relevant for the present study are, how- 
ever, the differences between converted and unconverted photons. 

According to the results of the preceding Section, we expect the opposite be- 
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atmospheric depth X (g cm" ) 

Fig. 4. Examples of shower profiles (randomly selected) for 10^^ eV photons at 
Auger North (solid black line), at Auger South (solid red line), protons (dashed 
green line), and iron nuclei (dotted blue line). For all events, arrival directions are 
fixed to local geographic North at 45° zenith. 

havior (larger Xmax at Auger North) for certain other directions. Results of 
simulations with the same conditions as in Figure H] and in the upper panel of 
Figure [5l but with the azimuth changed by 180°, are presented in the lower 
panel of Figure [51 The angle between this direction and the local magnetic 
field vector at Auger North is relatively small. As expected, the X^ax dis- 
tribution of photons at Auger South is now peaked at smaller values, while 
the distribution at Auger North is dominated by the large X^ax values from 
unconverted events. 

Another important quantity to distinguish photons from hadron primaries is 
the muon content of the air shower. It is well known that photon showers 
have less muons compared to hadron-initiated ones; important for the current 
investigations is whether a significant difference in the muon content exists 
between converted and unconverted photons. In Figure [6], for the same events 
shown in Figure HI the longitudinal profiles of muons are displayed. One can 
see a shift between the muon profiles for Auger North (all converted here) and 
Auger South (all unconverted here). This shift is, to a large extent, connected 
to the shift in Xmax (Figure H]). 

The distribution of the total number of muons at ground is given in Figure [71 
As expected, fluctuations are larger in case of unconverted events, i.e. at Auger 
South in this example. However, average values are quite similar. Thus, dif- 
ferences between converted and unconverted photon showers in the number 
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Fig. 5. Upper panel: Xmax distributions of different primaries for the same condi- 
tions as in Figured] (photons arriving from local geographic North). Lower panel: 
same as in upper panel, but the photon simulations were performed with the az- 
imuth changed by 180° (photons arriving from local geographic South). If indicated, 
distributions were scaled. 



of muons are much less spectacular than in X^ax- Part of the differences can 
directly be related to X^ax- We conclude that the main difference between 
converted and unconverted photons, and thus the main difference that could 
be exploited by observing from two sites, can be characterized by the Xmax of 
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Fig. 6. Examples of muon profiles (randomly selected) for lO^'^ eV primaries for a 
muon energy threshold 1 GeV (same events as shown in Figure U]) . 
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Fig. 7. Distributions of total number of muons at ground (energy threshold 1 GeV) 
for different primaries for the same conditions as in Figures S] and [S] (upper panel) 
and [H If indicated, distributions were scaled. 



the shower. 



There are other observables from ground arrays that were shown (or are ex- 
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pected) to differ between converted and unconverted photons, namely the 
shower size (total number of electrons at ground) [IT], the curvature of the 
shower front and the steepness of the lateral distribution of ground parti- 
cles [13], or the risetime of the detector signal at a certain distance from the 
shower core [?]. The differences between these observables for the two classes of 
photon events are closely connected to the stage of shower development and, 
correspondingly, to X^ax- We therefore expect, firstly, that also data from 
ground arrays alone, even without a direct observation of X^ax such as from 
fluorescence telescopes, can be exploited to search for complementary photon 
signatures at the two sites. And secondly, a study of X^ax distributions (which 
provides a measure of the relative shower age at ground) provides us with most 
relevant information for investigating possible complementary features of both 
sites. 



4 UHE photon scenarios and their observation at Auger South and 
North 

We now give an illustration of how advantage could be taken of the comple- 
mentarity between the preshower characteristics at Auger North and South. 
As examples, we consider a diffuse photon signal, a signal from a source region, 
and the absence of photons (i.e. determining upper limits). 



4.1 Diffuse photon signal 

In this example, we assume an isotropic primary flux with the all-particle 
energy power low spectrum with index -2.84. Such a spectrum is consistent 
with the first estimate from the Auger South Observatory [25]. It should be 
noted that a deviation from this power law, for instance a suppression of 
the flux above the energy of the speculated GZK cut-off, would result in a 
correspondingly modified all-particle flux at highest energies. The existence of 
such a hypothetic suppression will not affect the conclusions of this section: 
discussed differences between photon showers observed at Auger North and 
South depend only on the local geomagnetic field (see below for a quantitative 
example of applying a simple cut-off to the energy spectrum). 

We assume protons and photons as primaries. The input fraction of photons 
as a function of primary energy follows the results from a topological defect 
model in [2]. More specifically, the photon fraction is 9%, 15%, 27%, and 48% 
for energies above \g{E/eV) = 19.0, 19.3, 19.6, and 20.0, respectively (see also 
the curve labeled "TD" in Figure 8 in Ref. [7]). 
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For each Auger site, we simulated ~1000 events above lg(£'/eV) = 19.6 with 
zenith angles between 30—75° and random azimuth. We accounted for a de- 
tector resolution of 25 g cm~^ in X^ax and 10% in primary energy. The zenith 
angle range was chosen similar to that in the analysis of Auger data for the 
first photon limit [7j. In particular, near- vertical events are excluded in this 
example to reduce the fraction of photon events with Xmax below ground (see 
also discussion in Section [4.3p . The simulated event number corresponds to 
an observation time of about 5—8 years of the completed Auger South ground 
array. For a similar data taking period, data statistics of the fluorescence de- 
tectors is reduced by a factor ~7. 

In Figure [8] we show, as a scatter plot, the X^ax versus primary energy. We 
restricted the azimuth range in this plot to the local northern sky by requir- 
ing an azimuth between 30—150°. In this region of the sky, photon conversion 
starts at Auger North at smaller energies than at Auger South. As can be seen 
from the Figure, there are considerably less events with large X^ax (e.g. ex- 
ceeding 1000 g cm~^) at Auger North, for the same overall cuts applied to the 
data at both sites. 

The corresponding average X^ax and the X^ax fluctuations as a function of 
energy are shown in Figure [H Clear differences between the two sites can be 
noted. In particular, the larger average X^ax at Auger South is accompanied 
by significantly increased shower fluctuations. 

The Xmax distribution of events (as projected from the scatter plot shown in 
Figure [8]) is given for 50—100 EeV primary energies in Figure [TOl At Auger 
North, the events with large Xmax (unconverted photons, present at Auger 
South), are mostly lacking. Instead, events accumulate at ~ 900 g cm~^ (con- 
verted photons). To see the effect of a hypothetic GZK suppression of the 
energy spectrum at high energies, we applied a simple cut-off at E = 80 EeV 
(lg(-E/eV) = 19.9) to the data set. As expected, all the characteristic peaks 
of the Xmax distribution remain at the same positions as in case of spectrum 
without cut-off plotted in Figure [TOl The only difference is the heights of the 
peaks corresponding to converted and unconverted photons with respect to 
the proton peak: they decrease by 10-20% in case of spectrum with cut-off. 
This is due to the larger fraction of photons at higher energies in the assumed 
model. 

An observation of such different characteristics at Auger North and South 
(signal at different values of X^ax, behaviour of average X^ax and X^ax fluc- 
tuations as a function of energy), would be an unambiguous confirmation of 
a photon signal detection. 

So far, the local northern sky was studied. As expected, a similar but opposite 
behaviour is found for the local southern sky (not shown) when selecting, for 
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Fig. 8. Xjnax versus energy for events arriving from the local northern sky (see text 
for definition) at Auger North (top) and at Auger South (bottom), for unconverted 
photons (open blue circles), converted photons (filled black squares), and protons 
(red dots). 



instance, shower arrival directions from azimuth angles between 210—330°. 
Since this occurs at a somewhat higher energy, the event statistics is reduced 
in this case due to the steep primary flux spectrum. 
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Fig. 9. Average Xmax (top) and rms of X^az.^ (bottom) versus energy for the condi- 
tions shown in Figure [8l i.e. local northern sky at Auger North and at Auger South. 
For comparison, values corresponding to a pure proton flux are also shown using 
the model QGSJET 01 (dotted blue hues). 



Photon signal from a localized source region 



Let's now regard a hypothetical scenario of a photon signal coming from a 
certain source region. Localized hadron-accelerating sources are expected to 
produce photons at an energy typically 10 times smaller than the maximum 
energy of accelerated hadrons. Detection of the currently highest energy event 
at 3 X 10^° eV provides a motivation for the investigation of photons at not 
much lower energies produced by a localized source. 

The sky region accessible to both Auger sites is limited to source declinations 
<~ 53°. Requiring a maximum zenith angle of 75° for shower observations 
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Fig. 10. -'^max distributions for Auger North and South at the local northern sky 
at 50—100 EeV (projection from the scatter plot shown in Figure [8]). The different 
primary types (protons, converted and unconverted photons) contribute mainly to 
ranges in X^ax as indicated in the plot. In particular, the peak in the distributions 
at Xjnax ~ 800 g cm~^ comes from protons, at ~ 900 g cm~^ (at Auger North) 
from converted photons, and at ~ 1050 g cm~^ (at Auger South) from unconverted 
photons. 

at both sites, the range of accessible declinations is \6\ <~ 38°. 

The conversion probability of UHE photons from the direction of a certain 
source can change with time even at one site only, depending on the source 
declination. Examples are given in Figure [TT] for 10^^'^ eV and 10^" ° eV pho- 
tons. For instance, for a source at the celestial equator conversion probabilities 
between 0—75% are covered at 10^^'^ eV at Auger North. 

It is interesting to note that for the same conditions, the conversion probability 
is negligible at Auger South. Thus, characteristic differences when observing 
events from the same sky region exist that again could be used for a com- 
plementary study of UHE photons. A special case is the scenario of a source 
region observed at the same time from Auger North and South. Also here, 
differences exist in the instantaneous conversion probability of photons. This 
can be seen in Figure [TT] by comparing the curves at Auger North and South 
for the same source declination as a function of time. For instance, the instan- 
taneous difference in the conversion probabilities for 10^° '' eV photons from 
a source at the celestial equator can reach, in the extreme case, almost 100%. 
Differences exist also for sources at declinations S ~ —30°, i.e. for the galactic 
center region (Figure fTTil . 
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Fig. 11. Conversion probability at the two Auger sites versus time (in hours, offset 
arbitrary) for 10^^'^ eV (top) and lO^'^'^ eV (bottom) photons for different source 
dechnations as indicated. The same right ascension is assumed for all sources. The 
heavy lines denote the time intervals during which the sources are visible at zenith 
angles less than 75°. For better visibility the beginnings of lines are marked with 
large dots. 

4-3 Absence of photons 

If the photon flux is below the detection limits of current experiments, upper 
limits are deduced by comparing the data to expectations for showers initiated 
by primary photons. The signatures from photon showers differ between Auger 
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North and South due to geomagnetic cascading. Thus, even in case of observing 
no photons at both sites, differences arise in the analysis of the data since the 
photon hypothesis tested when deriving a hmit is different at each site. Upper 
hmits obtained at Auger South can, to some extend, be cross-checked this way 
at Auger North. 

In photon hmit studies, simulations are used to optimize the analysis cuts and 
to calculate the acceptance to photons. Due to the generally larger Xmax, the 
acceptance to photons can be smaller than that to nuclear primaries when 
requiring the X^ax to be above ground. Then, an efficiency correction to the 
photon limit may be needed. To keep the correction small, a fiducial volume cut 
was applied in [7] , requiring a minimum zenith angle of 35-42 deg (depending 
on energy). 

This condition could be relaxed at Auger North due to the onset of the 
preshower effect at smaller energies. On the one hand, the power of shower- 
to-shower discrimination is reduced for events from these arrival directions, 
as converted photon showers are more similar to proton-induced showers than 
unconverted ones. On the other hand, the result is a gain in instantaneous 
aperture at Auger North. The effective gain is largest at energies around 
lg(i?/eV)~19.8 where for a sizeable part of the sky, photons already can con- 
vert at Auger North, contrary to Auger South. For identical experiments and 
requiring Xmax to be above ground, the aperture is increased at Auger North 
by ~30% (20%) when allowing zenith angles up to 60° (75°) at lg(E/eV)~19.8. 



5 Conclusion 

The difference of conversion probabilities is a convenient and effective param- 
eter to estimate the complementarity between two sites for searching for UHE 
photons. Regarding the two sites of the Pierre Auger Observatory, Auger North 
and Auger South, significant differences in the preshower features of UHE pho- 
tons exist. The sky patterns of the preshower features are shifted due to the 
different pointing directions of the local magnetic field vectors. In addition. 
Auger North is located closer to the magnetic pole. This leads to a larger 
transition range in energy from small (< 10%) to large (> 90%) conversion 
probabilities compared to Auger South. More specifically, photon conversion 
starts at Auger North at smaller energies due to the larger magnetic field, 
but conversion probabilities > 90% are reached for the whole sky at higher 
energies due to the smaller bending of the field line with distance. 

These differences in the preshower characteristics result in different rates of 
(un-) converted photons from the same (both in local and astronomical coordi- 
nates) regions of the sky. Air showers initiated by converted and unconverted 
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photons can be well distinguished by current experiments. The main differ- 
ence is related to the position of depth of shower maximum Xmax? which is 
typically ~200— 300 g cm~^ smaller for converted photons. 

For a variety of UHE photon flux scenarios (diffuse photon flux; photons from 
source regions; absence of photons), the different preshower characteristics 
at the experimental sites can be used for a complementary search for UHE 
photons. Most important, a possible detection of UHE photons at Auger South 
may be confirmed in an unambiguous way at Auger North by observing the 
well predictable change in the signal from UHE photon showers. 

Acknowledgements: This work was partially supported by the Polish State 
Committee for Scientific Research under grants No. PBZ KBN 054/P03/2001 
and 2P03B 11024, by Polish Ministry of Science and Higher Education un- 
der grant N202 090 31/0623 and in Germany by the DAAD under grant 
No. PPP 323. Helpful remarks of the unknown referee, especially those con- 
cerning the motivation of Subsection 14. 2^ are kindly acknowledged. MR ac- 
knowledges support from the Alexander von Humboldt foundation. 



A Transverse component of the geomagnetic field along various 
trajectories 

B Preshower sky maps 

References 

[1] P. Bhattacharjee, G. Sigl, Phys. Rep. 327, 109 (2000); S. Sarkar, Acta 
Phys. Polon. B35, 351 (2004). 

[2] G. Gelmini, O.E. Kalashev, D.V. Semikoz, |arXiv:astro-p h/0506128]. 

[3] M. Ave et al, Phys. Rev. Lett. 85, 2244 (2000); Phys. Rev. D65, 063007 (2002). 

[4] K. Shinozaki et al, Astrophys. J. 571, L117 (2002). 

[5] M. Risse et al, Phys. Rev. Lett. 95, 171102 (2005). 

[6] G.I. Rubtsov et al, Phys. Rev. D 73, 063009 (2006). 

[7] Pierre Auger Collaboration, submitted to Astropart. Phys.; 
arXiv:astro-ph / 0606619 1 . 

[8] J. Abraham et al., P. Auger Collaboration, Nucl. Instrum. Meth. A 523, 50 
(2004). 



20 



0.1 



0.01 



0.1 



0.01 



AN, 0° 


^=90° (geograpliic Nortli) 

AN,'30°~ / - , AN, 60° 


AS, 0° 




: AS, 30° 
X AS, 6O7 




1000 10000 

altitude [km] 


<^=270° (geographic Soutli) 




~ ~ ~ 60° 
AN,60°~ ^ AS,30^\\ : 


AN, 30° 









1000 

altitude [km] 



10000 



Fig. A.l. Transverse component of the geomagnetic field (according to the IGRF 
model) along various directions at two observation sites: Auger South (AS) and 
Auger North (AN). The directions are specified by the zenith (attached to each 
curve) and azimuth (the same within one plot) angles given in the local coordinate 
system. The local minima occur at the altitudes at which the geomagnetic field is 
nearly parallel to the considered directions. 

[9] B. McBreen, C.J. Lambert, Phys. Rev. D 24, 2536 (1981). 

[10] F.A. Aharonian, B.L. Kanevsky, V.V. Vardanian, Astrophys. Space Sci. 167, 



21 




AN- AS 
10^3''eV 



















20% 














\ 80% 






')) 







,0% 



^0% 



w 



Fig. B.l. Sky maps of conversion probabilities at lg(£'/eV)=19.4 (left) and 19.6 
(right) for Auger North (AN, top), Auger South (AS, middle), and the difference 
(bottom); see also Figured] for further explanations. Additionally shown as dotted 
black lines are the declinations as projected on the sky for a given site (top and 
middle rows of plots) with a stepsize of 10°, with the celestial equator (declination 
= 0) highlighted as a dashed line and the pole indicated as a dot. 

Ill (1990); H.P. Vankov, P.V. Stavrev, Phys. Lett. B 266, 178 (1991); 
F. Halzen et al., Astropart. Phys. 3, 151 (1995); S. Karakula, W. Bednarek, 
Proc. 24^^ Intern. Cosmic Ray Conf., Rome, 266 (1995); A.V. Plyasheshnikov, 
F.A. Aharonian, J. Phys. G 238, 267 (2002); W. Bednarek, New Astronomy 
7, 471 (2002); H.P. Vankov, N. Inoue, K. Shinozaki, Phys. Rev. D 67, 043002 
(2003); P. Homola et al, Proc. 28*'" Intern. Cosmic Ray Conf., Tsukuba, 547 
(2003); H.P. Vankov et al., Proc. 28*^ Intern. Cosmic Ray Conf., Tsukuba, 523 
& 527 (2003); M. Ave et al, Phys. Rev. D 67, 043005 (2003); P. Homola et 
al., Acta Phys. Polon. B 35, 1893 (2004); M. Risse et al, Astropart. Phys. 21, 
479 (2004); A.C. O'Neill and S. Westerhoff, HiRes Collaboration Proc. 29*'' 



22 



AN, E AN, 10^°°eV 




Fig. B.2. Sky maps of conversion probabilities at lg(£'/eV)=19.8 (left) and 20.0 
(right) for Auger North (top), Auger South (middle), and the difference (bottom); 
see also Figures [T] and IB. II for further explanations. 

Intern. Cosmic R ay Conf., Pune 7, 373 (2005); M. Risse et al, Czech. J. Phys. 
56 (2006) A327, [ arXiv:astro-ph/0512 434 1 . 

[11] T. Stanev, H.P. Vankov, Phys. Rev. D 55, 1365 (1997). 

[12] W. Bednarek (1999), |arXiv:astro-ph/9911266] . 

[13] X. Bertou, P. Billoir, S. Dagoret-Campagne, Astropart. Phys. 14, 121 (2000). 
[14] P. Homola et al, Comp. Phys. Comm. 173, 71 (2005). 
[15] T. Erber, Rev. Mod. Phys. 38, 626 (1966). 

[16] T. Pierog et al, Nucl. Phys. B (Proc. Suppl.) 151, 159 (2006); T. Bergmann et 
al., [arXiv:astro-ph/0606564j . 



23 



AN, 10^°-^ eV E AN, lO^"" eV 




Fig. B.3. Sky maps of conversion probabilities at lg(£'/eV)=20.2 (left) and 20.4 
(right) for Auger North (top), Auger South (middle), and the difference (bottom); 
see also Figures [T] and IB. II for further explanations. 

[17] National Geophysical Data Center, USA, jhttp://www.ngdc. noaa.gov 

[18] N. A. Tsyganenko, National Space Science Data Center, 

NASA GSFC, Greenbelt, MD 20771, USA 

http: / / nssdc.gsfc.nasa.gov/space/model/ magnetos / data-based / geopack.html 

[19] T. Abu-Zayyad et al, HiRes Collaboration, Astropart. Phys. 23, 157 (2005). 

[20] M. Fukushima et al, Prog. Theor. Phys. Suppl. 151, 206 (2003). 

[21] L.D. Landau, I.Ya. Pomeranchuk, Dokl. Akad. Nauk SSSR 92, 535 & 735 
(1953); A.B. Migdal, Phys. Rev. 103, 1811 (1956). 

[22] D. Heck et al, Reports FZKA 6019 & 6097, Forschungszentrum Karlsruhe 
(1998). 



24 



[23] S. Eidelmann et al, Particle Data Group, Phys. Lett. B592, 1 (2004). 

[24] N.N. Kalmykov, S.S. Ostapchenko, A.I. Pavlov, Nucl. Phys. B (Proc. Suppl.) 
52B, 17 (1997). 

[25] P. Sommers for the Auger Collaboration, Proc. 29*^ Intern. Cosmic Ray Conf., 
Pune, 7, 387 (2005); [arXiv:astro-ph/0507150] . 



25 



